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PREFACE  AND  SUMMARY 

This  Memorandum  is  one  in  a  continuing  series  of  RAND 
•k 

publications  dealing  with  theoretical  and  computational 
questions  which  have  arisen  in  connection  with  the  RAND 
program  of  research  on  rocket  fuels  and  propellents  and 
in  biology  and  physiology.  It  details  a  method  by  which 
a  multi-  or  single-phase  chemical  equilibrium  problem  with 
a  large  number  of  different  chemical  species  can  be  re¬ 
placed  by  one  with  relatively  few  species.  The  purpose 
of  this  research  is  to  find  methods  for  transforming  com¬ 
plex  problems  into  problems  which  are  more  amenable  to 
numerical  solution. 

This  Memorandum  is  concerned  with  a  technique  first 
developed  in  an  earlier  Paper,  On  the  Reduction  of  Certain 
Multiplicative  Chemical  Equilibrium  Systems  to  Mathematical! 
Equivalent  Additive  Systems,  P-2419,  by  G.  B.  Dantzig  and 
J.  C.  DeHaven.  The  present  Memorandum  generalizes  the 
technique  to  systems  more  complex  than  those  to  which  the 
technique  of  the  earlier  paper  is  readily  applicable.  The 
basic  approach  of  the  previous  paper  is  described  briefly, 

Vc 

For  example,  see  Refs.  2-8. 


iv- 


but  the  reader  should  be  familiar  with  P-2419  for  certain 

explanatory  and  illustrative  material  which  has  not  been 
included  in  the  present  paper. 


-1- 


A  GENERALIZED  TECHNIQUE  FOR  ELIMINATING  SPSCIES  IN 
COMPLEX  EQUILIBRIUM  CALCULA.TIONS 

I.  INTRODUCTION 

We  will  consider  multi-  or  single-phase  chemical 
equilibrium  problems,  showing  that,  under  certain  circum¬ 
stances,  a  problem  involving  a  very  large  number  of  dif¬ 
ferent  chemical  species  may  be  replaced  by  one  involving 
relatively  few  species.  The  new  problem  will  then  be 
much  more  amenable  to  numerical  solution  than  the  original 
problem.  We  will  use  the  same  approach  in  formulating 
complex  chemical  equilibrium  problems  as  was  used  in  Cl]. 

It  is  briefly  described  here  primarily  in  order  to  intro¬ 
duce  our  notation  and  terminology. 

Consider  a  problem  involving  certain  chemical  species, 
X- , . . . ,X  , . . . ,X^ ,  in  one  or  more  phases.  It  is  no  loss 
of  generality  to  assume  that  no  species  can  occur  in  more 
than  one  phase.  For  example,  if  H2O  occurs  in  both  a 
liquid  phase  and  a  gaseous  phase,  we  assume,  as  is  fre¬ 
quently  done,  that  we  are  dealing  with  two  distinct  species, 
H20-liquid  and  H20-vapor,  and  we  formulate  the  mass  balance 
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equations  (see  below)  in  such  a  way  that  H^O-liquid  is 
transferable  into  or  from  H^O- vapor.  Thus,  we  assume  that 
species  in  distinct  phases  are  distinct. 

In  [l]  it  was  shown  that  any  chemical  equilibrium 
could  be  formulated  from  the  following  three  types  of 
information: 

I.  A  specification  of  which  species  occur  in  which 

phases . 

II.  Certain  constants  , . .  .  For  a 

discussion  of  the  meanings  and  the  various  interpretations 
of  these  constants,  as  well  as  the  methods  by  which  they 
may  be  obtained,  see  [2].  Because  one  way  of  looking  at 
the  depends  on  the  Gibbs  free  energy  function,  we 

will  call  the  free  energy  parameter  associated  with 

X^.  In  fact,  . 

III.  A  set  of  N  linear  equations  of  the  form 

L 

I  “  =0,n  .  1  ^  "  N  (1) 

t=l 

where  x.  is  the  (unknown)  number  of  moles  of  X..  Con- 
straints  of  the  type  (1)  may  represent,  for  example,  mass 
balance  constraints,  or  may,  for  example,  express  the 
condition  that  a  certain  phase  be  electrically  neutral 
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Thus,  if  n  of  Eq.  (1)  is  a  mass  balance  constraint, 

3.  might  be  the  number  of  atoms  of  some  chemical  element 
“V  ^  n 

in  each  molecule  of  X,,  and  s^  would  be  the  total  number 

I  0,n 

of  atoms  (measured  in  moles)  of  that  element  in  the  system 
under  consideration.  In  practice,  of  course,  it  is  ordin~ 
arily  more  convenient  to  base  the  mass  balance  equations 
on  various  radicals  rather  than  on  individual  atoms.  In 
any  case,  all  the  mass  balance  constraints  are  of  the 
form  (1);  that  is,  they  are  linear  equations  in  the  x  's. 

If  equation  n  of  (1)  is  an  electro-neutrality  con¬ 
dition  on  some  phase,  then  s.  would  be  zero,  except  for 

^  j  11 

those  I  for  which  X  occurs  in  that  phase,  and  X  is  an 

ion.  In  that  case,  s.  is  the  valence  of  X,,  with  the 

^ ,  n  -i 

appropriate  sign. 

If,  for  each  0  ^  ^  L,  we  define  a  vector,  S  , 

whose  components  are  s  ,  1  ^  n  ^  N,  then  Eq.  (1)  may  be 

n  j  v 

written  as 


L 


1=1 


(!') 


Notice  then,  that  is  essentially  the  "empirical  formula" 
for  X^. 

We  will  call  the  constraint  vector  associated  with  X^ . 
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Repeating,  a  chemical  equilibrium  problem  may  be 
completely  specified  by  the  following  three  types  of  in¬ 
formation: 

I,  The  specification  of  which  species  occur  in 
which  phases, 

II.  An  energy  parameter,  for  each  species, 

III.  A  constraint  vector,  S^,  for  each  species, 
and  a  vector  Sq. 

Note  that  the  concentration  [X.  1  of  each  species 

0 

X  can  be  calculated  by  dividing  x  by  the  sum  of  all 
0  0 
X . ' s  for  which  X  occurs  in  the  same  phase  as  X  (this 

i  -t  Iq 

sum,  of  course,  includes  x.  itself).  Thus,  the  [X . ] ' s 
are  functions  of  the  * s  (although  each  [X^ ]  does  not 
depend  only  on  the  corresponding  x^) . 

It  was  shown  in  [2]  and  [3]  that  the  x  's  may  be 
efficiently  determined,  numerically,  by  minimizing  the 
function 

L 

1=1 


•k 

Throughout  this  paper,  we  measure  concentration  in 
mole  fractions. 
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subject  to  the  constraints  that 


L 


1  =  Sq  . 


(1’) 


Before  stating  and  proving  our  results,  we  give  an 
example  of  a  typical  and  relatively  si-mple  situation  to 
which  they  apply. 

Consider  a  chemical  equilibrium  problem  involving  one 
or  more  phases.  Consider  some  particular  liquid  phase, 
which  we  will  arbitrarily  call  the  first  phase.  Suppose 
that  a  certain  chemical  species,  X,  occurs  in  the  first 
phase.  Suppose  also  that  X  can  combine  with  oxygen  to 
form  new  species  XO2  and  XO^.  Suppose  that  X  has  50 
distinct  sites  at  which  hydrogen  ionization  can  occur, 
that  XO2  has  75  distinct  sites  at  which  hydrogen  ioniza¬ 
tion  can  occur,  and  that  XO^  has  100  distinct  sites  at 
which  hydrogen  ionization  can  occur. 

Suppose  that  all  these  sites  are  independent,  in  the 
same  sense  that  ionization  at  one  site  does  not  influence 
the  ease  or  difficulty  of  ionization  at  any  other  site 
on  the  same  molecule.  (Note:  This  assumption  of  complete 
independence  is  not  necessary  in  order  to  apply  the  methods 
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of  this  paper.  The>  ^ould  also  be  applicable  if  various 
types  of  systematic  dependence  were  assumed.) 

Suppose  also  that  all  these  sites  are  of  two  distinct 
classes--sites  of  the  first  class  being  easier  to  ionize 
than  sites  of  the  second  class.  But,  suppose  that  all 
sites  of  the  first  class  are  equally  easy  to  ionize  and 
that  all  sites  of  the  second  class  are  equally  easy  to 
ionize.  Thus,  each  molecule  of  X  might  have  20  sites  of 
the  first  class  and  30  sites  of  the  second  class;  each 
X0.J  molecule  might  have  30  sites  of  the  first  class  and 
45  sites  of  the  second  class.  Each  XO^  molecule  might 
have  50  sites  of  the  first  class  and  50  sites  of  the 
second  class. 

Notice  that  two  molecules  ionized  in  different  ways 

represent  different  species;  hence,  we  are  dealing  with 

a  total  of  2^^  different  varieties  of  X,  2^^  different 

varieties  of  XO2J  and  2^^^  different  varieties  of  XO^-- 

30 

or  a  total  of  approximately  10  distinct  chemical  species. 

It  will  turn  out  that  we  will  be  able  to  formulate  a  new 

problem  which  will  be  equivalent  to  the  first  problem, 

30 

but  one  in  which  all  these  approximately  10  species  will 
'e  replaced  by  only  seven  species.  Of  course,  the  first 
phase  may  itself  contain  additional  species  and  there  may 
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be  other  phases.  But  we  will  have  attained  a  significant 
reduction  in  problem  complexity. 

We  V7ill  now  leave  this  special  example  and  proceed 
to  state  and  prove  our  general  results.  Although  we  will 
not  again  refer  to  this  example,  the  reader  may  wish  to 
refer  back  to  it  from  time  to  time  in  order  to  understand 
the  applicability  of  the  general  and  rather  abstract  dis¬ 
cussion  which  follows.  It  should  be  emphasized  that  the 
example  is  not  representative  of  either  the  complexity  or 
the  variety  of  the  situations  to  which  our  results  may  be 
applied.  The  only  purpose  of  the  e;cample  is  to  enable 
the  reader  to  more  readily  understand  the  statement  of 


our  results. 
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II.  RESULTS 

Let  (])  be  some  subset  (or  the  entire  set)  of  species 
in  one  of  the  phases  (which  we  will  arbitrarily  call  the 
first  phase)  of  a  multi-  or  single-phase  equilibrium 
problem  in  the  sense  of  [l].  Let  (j)  be  partitioned  into 
various  sets  of  species,  .  ,0^.  That  is,  every 

species  9  of  (J  is  in  exactly  one  of  the  sets,  (jl^. 

Let  j .  .  .  ,A^, . . .  ,Aj^  be  classes  of  sites  (these  sites 
are  actually  abstract  entities,  and  need  not--but  may-- 
correspond  to  actual  chemical  sites) .  Suppose  that  a  site 
in  the  class  can  have  possible  states: 

B,  .  (For  example,  a  class  of  sites  might  be 

the  class  of  all  oxylabile  sites  which  are  affected  in  a 
given  manner  by  the  presence  of  O2;  the  possible  states 
of  a  site  in  this  class  might  then  be  an  ionized  state 
and  an  un-ionized  state.) 

Suppose  further  that  for  each  i  and  k,  every  species, 
9,  in  the  set  d)^  has  exactly  the  same  number  of  sites, 

Pik,  in  the  class  A^.  That  is,  the  classification  of  the 
sites  of  a  species,  9,  depends  only  on  which  of  the  sets, 
(j)^,  that  9  belongs  to. 
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Suppose  also,  that  any  species,  0,  of  (j)  is  completely 
identified  by  stating  the  of  which  it  is  a  member  and 
the  states  of  each  of  its  sites. 


Suppose  also,  that  for  each  0^,  all  logically  possible 
assignments  of  states  to  sites  (within  the  limitation  that 
if  a  site  is  in  the  class  then  its  state  must  be  one  of 
B, ,  , . . .  ,B,  . , . . . ,B,  -  )  actually  occur  among  the  species  of 

tCi  Kj 


Suppose,  finally,  that  the  free  energy  parameter,  ag, 
and  mass  constraint  vectors,  Sg,  for  each  0  of  (|>  depend 
only  on  which  of  the  (|)^'s  that  0  belongs  to  and  on  the 
states  of  the  sites  of  0;  and,  that  this  dependence  is  an 
additive  dependence.  To  state  this  assumption  more  pre¬ 
cisely,  for  each  species  0  of  ({),  and  for  each  site  class 
A^,  and  for  each  j,  1  ^  j  ^  J^,  let  h(0,k,j)  be  the  number 
of  sites  of  0  of  class  which  are  in  the  state  B^^ .  Our 
assumption  then  states  that  for  each  set  0^^  of  species 

there  is  a  vector  T.  and  a  number  i3.  ,  and  that  for  each 

1  1 

site  class  A^  and  each  j  ,  1  -  j  ^  there  is  a  vector 

T,  .  and  a  number  .  such  that 
kj  kj 


^0  =  \ 


K  ^ 

^  I  h(0,k,j)T^^  ,  for  all  0€0. 


k=l  j=l 


(2) 


-lo¬ 


an  d 


“9 


K  '"k 

^.  +  V  ^  h(9.k,j)^^  ,  for  all  9€((.  . 

k=l  j=l 


(3) 


We  may  now  state  our  results. 

Theorem  1.  Let  Xg  be  the  number  of  moles  of  species 
0  present  in  the  first  phase  at  equilibrium.  For  each  set 
0^  of  species  let 


L  U  I 


9 


(4) 


9C0. 


and  for  each  state  B,  .  let 


'kj 


=  V 


^  h(0,k,j)xg 
0€0 


(5) 


(Note  that  is  then  the  total  number  of  moles  of  species 

in  0.,  and  that  y,  .  is  the  total  number  of  sites  in  the 
^1’  "^kj 

state  3^^ --measured  in  moles-^regardless  of  the  species  on 
which  these  sites  occur.) 

Then,  the  following  relationships  hold: 


I: 

i=l 


1=^9  ’ 


9^0 


(6) 
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I 

1  Pik^i 

i=l 


for  all  site  cl/'sses,  A^; 


h(0,k,j) 


(7) 


(8) 


for  all  sets  and  all  species 


where  denotes 


Theorem  2.  If  ve  replace  the  original  problem  by  a 
new  problem,  as  follows: 

a)  Eliminate  all  the  species  0  of  from  the  first 
phase; 

b)  For  each  ().  ,  add  to  the  first  phase  a  new  species 

i  with  constraint  vector  and  energy  parameter 

B  Let  V  then  be  the  number  of  moles  of  species 
i*  i 

i  in  the  new  problem; 

c)  For  each  k,  1  ^  k  ^  K,  introduce  a  new  phase  whose 
species  will  be  the  states  ,  the  constraint 
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vector  and  energy  parameter  associaced  with  B,  . 
being  and  j  >  respectively.  Let  y^^  then 
be  the  number  of  moles  of  species  in  the  new 
problem; 

d)  All  the  species  of  the  original  problem  (except 
those  in  ())  will  be  retained  in  the  new  problem 
with  their  original  constraint  vectors  and  energy 
parameters ; 

e)  For  each  A^,  a  new  constraint,  Eq.  (7),  is  intro¬ 
duced  ; 

then,  the  two  problems  will  be  equivalent  in  the  following 
senses: 


I.  The  minimum  free  energy  of  the  original  problem 
is  the  same  as  the  minimum  free  energy  for  the  new  problem. 

II.  The  number  of  moles  of  all  the  species  (and 
their  concentrations),  except  for  those  species  in  (|l 
which  do  not  occur  in  the  new  problem,  will  be  the  same 
at  equilibrium  in  the  two  problems. 

III.  If  Xg  is  the  amount  of  species  0  in  the  first 
problem  at  equilibrium,  and  if  y^  and  y^^  are  the  amounts 
of  species  i  and  B,  .,  respectively,  in  the  new  problem  at 
equilibrium,  then  these  quantities  arc  related  by  Eqs. 
(4)-(S).  (Equation  (7)  holds,  of  course,  simply  because 
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it  was  imposed  as  a  constraint  on  the  new  problem.) 

Note:  It  should  be  observed  that  having  replaced  the 

old  problem  by  the  new  problem  and  having  solved  the  new 
problem,  we  may  recover  the  original  Xg's  by  direct  sub¬ 
stitution  in  Eq.  (8). 
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III.  PROOFS 


Throughout  this  section  we  will  assume  that  each 
and  each  is  a  given  vector  and  that  the  Sg's  are  given 
by  Eq.  (2),  We  will  also  assume  that  each  and  each 

is  a  given  number  and  that  the  QIq's  are  given  by  Eq.  (3) 


We  will  allow  the  y^'s,  y^^'s,  and  Xg’s  to  be  vari¬ 


ables.  We  will  make  explicit  various  different  assumptions 
about  the  manner  in  which  they  are  related  to  each  other 
in  the  various  lemmas.  Thus,  at  times  we  will  assume  that 
the  y’s  are  defined  in  terms  of  the  x's  by  Eqs,  (4)  and 
(5)j  and,  at  times  we  will  assume  that  the  x's  are  defined 
in  terms  of  the  y's  by  Eq.  (8). 

However,  we  will  always  regard  as  being  defined  by 


I 


k  L  \j  • 

j=l 

Lemma  1.  ^  h(0,k,j)  =  p^^,  for  all  0  in  the  set  (|)^ 

j=l 


and  for  all  k. 

Proof.  This  is  simply  a  restatement  of  the  assumption  that 
a  species  in  (j)^  has  p^^^  sites  in  the  class  of  sites. 

We  will  find  it  convenient  to  introduce  a  polynomial 
in  several  variables,  R(z) ,  where  the  components  of  z  are 
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z.,  for  and  z,  . ,  for  k=l,...,K  and  ; 

1  K.J  K 


K 


thus,  z  has  a  total  of  I  +  ^ components.  We  define: 


k=l 


K  '^k 


««■  z*.  z  n  n 


i=l  9€0^  k=l  j=l 


(9) 


Lemma  2. 


K  ‘^k 


,*-•  Z  M..M)  n  TT  ■ 


kj  ,1 


i=l  Oc0. 


k«l  j=l 


kj 


Proof.  By  calculation  from  9. 
Lemma  3. 


I  K  /  •’k  \ 

n  I  ^k j  • 

i=l  k=l  ^j=l  / 


Proof.  Note  that  we  have  assumed  that  all  the  0's  in 
the  same  0^  have  the  same  classification  of  sites.  Imagine 
that  for  each  0  of  0^  and  for  each  site  of  0  we  write 
(with  a  very  small  pen  I)  the  sum 
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if  the  site  is  in  the  class  A^.  Note  that  we  will  have 
written  each  sum  times  on  each  0  of  Let  us  then 

develop  the  product 


n 


(10) 


as  follows. 

For  each  0  of  write  a  product  of  z's  over  the 

sites  of  0  by  choosing  at  each  site  that  term  cor¬ 
responding  to  the  state  that  the  site  is  in.  Let  us 
then  add  the  resulting  products  for  all  0  of  On  the 

one  hand  it  is  clear,  from  the  definition  of  h(0,k,j),  that 
the  resulting  sum  will  be: 


K 


i:  n 


(11) 


k=l  j=l. 


On  the  other  hand,  since  we  have  assumed  that  as  0  varies 
over  (j>^  all  possible  permissible  assignments  of  sites  to 
states  occur  exactly  once,  it  is  not  bard  to  see  that  the 
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process  has  also  evaluated  the  product,  (10).  Thus,  the 
quantities  (10)  and  (11)  are  equal. 

If  (10)  and  (11)  are  set  equal  and  multiplied  by  z^, 
and  the  result  is  added  over  i,  we  arrive  at  the  identity 
asserted  in  Lemma  3. 

Lemma  4.  If  the  relationships  of  Eqs.  (4)  and  (5) 
hold,  then: 

I  K 

Z  + 

i=l  k-1  j=l  0c(ll 

Proof.  This  follows  by  expanding  Sg,  using  Eq.  (2),  and 
from  Eqs.  (4)  and  (5). 

Lemma  5.  If  the  relationships  of  Eqs.  (4)  and  (5) 
hold,  then 

I  K  ‘^k 

i  +  Z  Z  ykj\j  =  Z  ='0“0  • 

-1  k=l  j=l  0c(ll 

Proof.  From  (3),  (4),  and  (5). 

Recall  that  as  0  varies  over  (jl^,  each  possible 
assignment  of  states  to  sites  is  attained  exactly  once. 


Z  Z  =  Z  -0S, 
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Lemma  6.  If  Eqs.  (4)  and  (5)  hold,  then  so  does 
Eq.  (7). 

Proof.  Lemma  1. 

Lemma  7 .  If  Eq,  (4)  holds,  then  so  does  Eq.  (6). 

Lemma  8 .  If  the  x^'s  are  the  quantities  associated 
with  the  original  problem  and  if  the  y's  are  defined  by 
Eqs.  (4)  and  (5),  then  Eq.  (8)  holds. 

Proof.  We  know  that  for  an  appropriate  Lagrangian  vector, 

o(g  +  log  ~  • 

Here,  a  denotes  the  total  number  of  moles  at  equilibrium 
in  compartment  one  of  the  first  problem. 

Hence,  using  Eqs.  (2)  and  (3)  (and  letting  0  be  in  0^): 

X  K  ‘^k 

log  ^  =  „.Sg  -  ttg  =  (ir-T.-fi.)  +  I  I  h(9.k,j)(’r-\-^kj) 

k=l  j=l 


or,  letting 


=  e  ^  ^  >  0 
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and  letting 


>  0 


we  get 


K  '^k 


Xq  av^ 


i  n  n  .  for  eciti, . 


k=l  j=l 


kj 


(12) 


Let  c^,  1  ^  k  ^  K,  be  arbitrary  positive  numbers  and 


set  w,  . 
kj 


c,  V,  . . 

k  kj 


Also,  set 


w.  =  a 
1 


K 


k=l 


then,  using  Lemma  1,  Eq.  (12)  becomes 


X 


0 


K  "k 

1  r  ]  r  h(9,k,j) 

I  I  I  Kkj 

k=i  j=i 


for  0C{^. 

L 


Now,  by  taking 


(13) 
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1 


-21- 


=  (From  the  definition  of  R.) 


w. 

1 


Sw. 

1 


=  (From  (16))  . 


(18) 


Next,  by  using  Eqs.  (5)  and  (13),  we  get 


I 

E  "i  E 

i-1 


K 

k=l  j=l 


=  (By  Lemma  2)  w 


^3 


6R(w) 

6w,  . 


(19) 


Combining  (17) ,  (18) ,  and  Lemma  6 


6w,  . 


(20) 


Combining  (20)  and  (19) 


(21) 


then,  (18),  (21),  and  (13)  yield  (8). 
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Lemma  9.  If  Eqs.  (7)  and  (8)  hold,  so  do  ilqs.  (4) 
and  (5) . 

Proof.  Let  w,  ,  =  — — ^  and  let  w.  =  y..  Then 
-  kj  1 


Z  “kj  =  1 


(22) 


and  Eq.  (8)  becomes 


K  “^k 


]~~|  for  all  0c(|)^  .  (23) 

k=l  j=l 


Next 


K  •^k 


^h(0,k,j)xg=  ^  ^h(0,k,j)  I 


Q€(^. 


i=l  0€(|)^  k=l  3=1 


(By  Lemma  2)  w^.  . 

kj 


(24) 


By  Lemma  3  and  Eq.  (22) 


|RM=  y.  „  =  Vp  y 

9w,  .  L  ^ik  1  L  ^ik^i 


=  (By  Lemma  6)  Yi  .  = 
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This  and  Eq.  (24)  yield  (5). 
Finally,  from  (23) , 


K 


z-.--.  I  n  n 


0€(J^  0C0^  k=l  j=l 


=  (From  the  definition  of  R)  w,  t ^ 

1  ciw. 

1 


But,  Lemma  3  and  Eq.  (22)  yield  =  1.  This  and 

i 

the  fact  that  =  y^  yield  (4). 

Lemma  10.  (A)  If  Eqs.  (4)  and  (5)  hold  and  if  the 

Xg  are  those  associated  with  the  equilibrium  solution  in 
the  original  problem;  or  (B)  if  Eqs.  (7)  and  (8)  hold. 


then; 


K  "^k 


I 


I  I  Wi 

L=1  k=l  j=l 


(25) 


I 

A  ' 


9  + 


)■  i-.(^ 

'  •  -I 


.  +  log  — 
1  a 


K  ‘^k 


I  Z 


(26) 


k=l  j=i 
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where  by  "a''  we  mean  the  sum  of  all  amounts  of  all  species 
in  the  first  compartment  of  the  original  problem. 

Note  1:  By  Lemmas  6  and  8,  (A)  implies  (B) .  Hence, 
we  need  only  prove  that  (B)  implies  (25)  and  (26).  But, 
by  Lemma  9,  (B)  implies  Eqs.  (4)  and  (5)  and  by  Leiima  7 
Eqs.  (4)  and  (5)  imply  Eq.  (6).  Hence,  in  proving  Lemma 
10  we  may  assume  Eqs.  (4),  (5),  (6),  (7),  and  (8). 

Note  2:  By  Eq.  (6),  a  is  also  the  sum  of  all  the 
amounts  of  all  the  species  in  the  first  phase  in  the  new 
problem. 

Equation  (25)  follows  at  once  from  Lemma  4. 

To  prove  Eq.  (26),  observe  that  by  employing  Lemma  5 
and  Eq.  (26)  we  need  only  prove  that 

I 

1  ^9  ’‘q  =  Z 

9€(})  i=l 


E  I  ^kj  ^  • 

k=l  j=l 


(27) 


To  prove  Eq.  (27),  note  that  from  Eq.  (8)  we  get 


K 

log  Xg  =  log  ^ 

k=l 


^  h(0,k, j) 

j=l 
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for  0  in  0^,  and  hence 


I 

0c(j) 


Xq  log  Xq  = 


(  1 1°8  yi  1  ^9 

\i=l  0c(l». 


K  ‘^k 

-  I  I 

k=l  j=l 


log 


^  Ih(e.k,j)xgj 

K  ^  -  a.  / 


0C(1) 


K  ‘^k 


=  I  Yi  log  Yi  +  X  I  y^j  log  - 


i=l 


k=l  j=l 


which  is  Eq.  (27).  QED. 

Theorem  1  follows  from  Lemmas  7,  6,  and  8. 

To  prove  Theorem  2,  observe  that,  given  a  solution 
of  the  old  problem,  we  may  define  the  y's  by  means  of 
Eqs.  (4)  and  (5)  to  obtain  a  solution  to  the  new  problem. 
Conversely,  by  means  of  Eq.  (8)  we  may  obtain  a  solution 
to  the  old  problem  given  a  solution  to  the  new  problem. 
Lemmas  4,  6,  8,  9,  and  10  assure  the  validity  of  this 
process  and  of  the  various  assertions  of  Theorem  2. 
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